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Abstract 

Three independent predictions follow from postulating the existence of protons 
co-accelerated with electrons in extragalactic jets (i) multi-TeV gamma ray emission 
from nearby blazars, (ii) extragalactic cosmic ray protons up to ~ 10 20 eV, and (iii) 
extragalactic neutrinos up to ~ 5 x 10 18 eV. Recent gamma ray observations of 
Mrk 421 and Mrk 501 employing the air-Cerenkov technique are consistent with the 
predicted gamma ray spectrum, if one corrects for pair attenuation on the infrared 
background. Prediction (ii) is consistent with cosmic ray data, if one requires that 
jets are responsible for a at least a sizable fraction of the extragalactic gamma ray 
background. With kubic kilometer neutrino telescopes, it will be possible to test (iii), 
although the muon event rates are rather low. Neutrino oscillations can increase the 
event rate by inducing tau-cascades removing the so-called Earth shadowing effect. 



1 Introduction 

In the early days after the discovery of extragalactic radio sources, it was a 
widely held belief that the relativistic electrons responsible for the observed 
synchrotron emission are secondary electrons from pp-interactions of acceler- 
ated protons on ambient gas (1; 2). Protons and ions were known to be the by 
far dominant species in the observed cosmic rays, and this was expected to be 
mirrored at their acceleration sites. Nevertheless, the picture was soon given 
up, since it would require enormous amounts of target matter in the jets which 
is inconsistent with plausible energetics and the non-observation of emission 
lines and bremsstrahlung. An electron-positron composition of the jets was 
also suggested from theory claiming that the plasma feeding radio jets should 
be due to pair production in the ergosphere of a maximally rotating black hole 
(3). The dielectric properties of radio jets inferred from observations of weak 
Faraday rotation and circular polarization in the jet of 3C279 seem to lend 
further support to a light composition (4). However, this picture also bears 
several fallacies. 
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Target matter does not have to be present to obtain efficient cooling of ac- 
celerated protons. Photo-production of secondary particles and synchrotron 
emission can become important, if the proton energy is high enough, since the 
cooling rate for these processes decreases with energy. In fact, in a statistical 
acceleration process in which there is enough time and space to balance accel- 
eration energy gains against energy losses, it is an inevitable consequence that 
the protons reach ultrahigh energies (5; 6). Even if the original composition of 
the jet plasma were light, polluting baryons from the ambient medium would 
quickly take over most of the jet's momentum, so that the acceleration mech- 
anism must eventually tap baryonic kinetic energy. In fact, for gamma ray 
bursts, this is believed to be the crucial explanation for the fact that the burst 
energy is tapped only by shock fronts far away from the site of the original 
pair fireball (8). Therefore, it is the most natural assumption that the observed 
nonthermal emission is partly due to the accelerated electrons, and partly due 
to the accelerated protons (7). Inferences of the plasma composition based on 
measurements of polarization and Faraday rotation are based on the assump- 
tion that the electron distribution observed in the optically thin synchrotron 
regime traces down to lower energies which is known to lead to theoretical 
inconsistencies (9). Since the radiative properties are determined by the par- 
ticles with the highest energies but the dielectric properties by those with the 
lowest energies, there is a general mismatch in the conclusions drawn from 
observations sensitive to either regime and one must be careful with claims 
about jet composition. 

It may be a good advice not be too narrow-minded when confronted with 
a new observational result such as multi-TeV emission from blazars and to 
include as many independent facts as possible in its interpretation. If the 
gamma rays were indeed primarily of a hadronic origin, there are a number of 
corollaries which allow to falsify the claim, whereas the arguments for a purely 
leptonic origin of the gamma rays boil down to some version of Ockham's razor 
( "it is more economic to use the observed electrons to model the gamma ray 
emission"). I wonder whether this is enough to get around the symptomatic 
fact that there was no prediction of multi-TeV emission from leptonic models 
prior to the observations. Moreover, there are problems with leptonic models, 
such as the missing intrinsic curvature in the multi-TeV spectrum of Mrk 501 
and the surprisingly low magnetic field values as was pointed out in ref. (10). 
Although the most recent HEGRA spectra of Mrk 501 (11) do show some 
curvature, the lower limits on the infrared background imply that the intrinsic 
spectrum must be rather flat or even up-turning (12). 

One might as well argue that it is much more in the sense of Ockham's ra- 
zor (i.e. more economic), if one finds that the same sources which very likely 
produce most of the extragalactic gamma rays would at the same time pro- 
duce cosmic rays at ultrahigh energies where it is difficult to find any other 
astrophysical source supplying enough power. There are two independent ob- 
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servations and only one model. According to Landau, a third independent 
fact is needed to consider a theoretical claim seriously. Here I consider the 
high-energy neutrino emission intimately connected with multi-TeV gamma 
rays from hadronic accelerators as the missing piece of information. To first 
elucidate the connection between extragalactic high-energy emissions in the 
framework of global energetics, Sect. 2 quantifies the non-thermal energy that 
may be released by active galactic nuclei and their jets integrated over their 
cosmic history. This qualifies extragalactic jets as possible sources of the ul- 
trahigh energy cosmic rays and implies multi-TeV emission from the jets due 
to proton energy losses at their acceleration site. Section 3 discusses the pre- 
dicted spectra from a simple quasi-stationary unsaturated synchrotron cascade 
emission model. Bearing on the assumption of a strong evolution of their lu- 
minosity density, the neutrino and cosmic ray spectra from extragalactic jets 
using the assumptions of the hadronic gamma ray emission model are com- 
puted in Sect. 4. Since the expected peak of the neutrino spectrum lies at 
energies in excess of 100 TeV at which the Earth becomes optically thick with 
respect to neutrino absorption, Sect. 5 extends the discussion by including the 
profound effects of neutrino oscillations. 



2 Origins of extragalactic background radiation 



Inspection of Fig. 1 shows an interesting pattern in the present-day energy 
density of the diffuse isotropic background radiation consisting of a sequence 
of bumps each with a strength that is decreasing with photon energy. The 
microwave bump is recognized as the signature of the big bang at the time 
of decoupling with its energy density given by the Stefan-Boltzmann law 
U3K = cT A . The bump in the far- infrared is due to star formation in early 
galaxies, since part of the stellar light, which appears as the bump at visible 
wavelengths, is reprocessed by dust obscuring the star-forming regions. The 
energy density of the two bumps can be related to the present-day heavy el- 
ement abundances. Heavy elements with present-day mass fraction Z = 0.03 
were produced in early bursts of star formation by nucleosynthesis with ra- 
diative efficiency e = 0.007 yielding the present-day energy density 

P*Zec 2 

1 + Zf 



where p* denotes the mass density of baryonic matter and Zf the formation red- 
shift corresponding to the era of maximum star formation. This is, of course, 
only a very rough approximation of the true star formation history, but good 
enough to set the scale for an argument pertaining to the global energetics. 
In particular, the ratio between the energy released by stars and by other 
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Extragalactic radiation background 




Fig. 1. Sketch of the present-day energy density of the extragalactic radiation back- 
ground from radio waves to gamma rays. 

sources with the same formation history is independent of its details. Let 
denote the baryon density in terms of the critical density of the Universe 
and h = H o /100 km s _1 Mpc -1 the dimensionless Hubble constant, then the 
energy density takes the value 

^ ~ 6 x 1CT 3 (^Q (i±*) " eV cm" 3 (2) 



and should represent the energy density of the sum of the far-infrared and 
optical bumps. Probably all galaxies (except dwarfs) contain supermassive 
black holes in their centers which are actively accreting over a fraction of 
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released by the accreting black holes amounts to 



^agnA* ~ 10 of their lifetime implying that the electromagnetic radiation 



eaccr^bh ^agn ^ _ 1 , m -j A - / 

ZeM* U 



-u ns ~ 1.4 x 10" 4 eV cm" 3 (3) 



adopting the accretion efficiency e accr = 0.1 and the black hole mass fraction 
Mbh/M* = 0.005 (13). Most of the accretion power emerges in the ultraviolet 
where the diffuse background is unobservable owing to photoelectric absorp- 
tion by the neutral component of the interstellar medium. However, a fraction 
of M x /M bh ~ 20% taken from the average quasar spectral energy distribution 
(14) shows up in hard X-rays producing the diffuse isotropic X-ray background 
bump with u x ~ 2.8 x 10 -5 eV cm -3 (15). Jets with non-thermal 7-ray emis- 
sion show up only in the radio-loud fraction £ r i ~ 20% of all AGN and their 
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kinetic power roughly equals the accretion power (16). Hence one obtains for 
the energy density due to extragalactic jets 

u i = (h) u ^ ~ (h) 2 - 8 x 10-5 eV cm " 3 - ( 4 ) 



This energy is released in relativistic particles, magnetic fields, and pdV ther- 
modynamic work against the ambient medium into which the jets propagate. 
Adopting a radiative efficiency of £ ra d = 10% for the jets, the gamma ray 
(cosmic ray, neutrino) energy released by the jets can amount to a maximum 
present-day energy density of 

u, ~ ~ 2.8 x 10- 6 eV cm" 3 (5) 



which comes remarkably close to the energy density 3.2 x 10 -6 eV cm -3 of the 
extragalactic gamma ray background observed between 100 MeV and 30 GeVQ 
using the spectrum given in ref. (17). Note that this is consistent with having 
a particle acceleration efficiency in radio jets which is of the same order of 
magnitude as the 13% efficiency required for supernova remnants to produce 
the Galactic cosmic rays. 

Protons can achieve a high radiative efficiency due to photo-production of 
secondary particles or synchrotron emission only if they reach ultrahigh en- 
ergies. This can be seen most easily for synchrotron radiation in which case 
the proton Lorentz factor must be larger than the electron Lorentz factor by 
(m p /m e ) 3 to produce the same cooling rate. Particles with energies exceed- 
ing 6 x 10 18 eV may seem outrageously exotic to many astrophysicists, but 
they are actually observed in the local cosmic ray spectrum. These particles 
cannot be confined by the Milky Way, but nevertheless show a near-isotropic 
sky distribution. Furthermore, the local cosmic ray spectrum flattens above 
6 x 10 18 eV indicating a separate extragalactic source population. It will be 
argued in Sect. 4 that only an extragalactic source population as strong as the 
one supplying the diffuse isotropic gamma rays can be considered for their ori- 
gin. Before highlighting these arguments, it is shown in the next section that 
protons at such high energies produce interesting gamma ray spectra owing 
to the photo-production of secondaries. 



Note that the flux in the gamma ray background observed by CGRO is close to 
the bolometric gamma ray flux, since pair attenuation and cascading must lead to 
a turnover of the background spectrum above 20 — 50 GeV for extragalactic source 
populations (18; 19) 
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3 Proton blazar predictions and observed multi-TeV spectra 

A large number of emission models for radio jets based on a shock-in-jet sce- 
nario exist and have been shown to explain most of the low-frequency jet 
phenomenology using a number of parameters such as jet kinetic energy, mag- 
netic field gradient, opening angle of the flow channel, etc. within plausible 
ranges. It is straightforward to include relativistic protons in such models as- 
suming they have a spectrum with the same slope but different maximum and 
minimum energies (due to different energy losses, gyro-resonant thresholds, 
and Larmor radii). Among the interesting consequences of the protons are 
(i) a higher nonthermal pressure and (ii) more gamma ray flux (adding to the 
Compton flux from the accelerated electrons). The equipartition magnetic field 
strength increases according to B oc k 2 ^ 7 where k = u p /u e denotes the ratio 
between the energy density in relativistic protons and electrons, respectively, 
and this may help to alleviate some problems for pair jets such as the observed 
pressure support of the radio lobes in NGC 1275 against the surrounding hot 
intracluster medium of the Perseus cluster (20). 

In terms of simplicity and predictive power, the original Blandford & Konigl 
conical jet model (21) is useful and has been investigated for the radiative 
signatures of the protons (7). The model describes the stationary emission from 
a conical section of a free relativistic jet and is therefore of limited applicability 
to non- stationary features in the observed spectra. The relativistic particle 
and magnetic energy density decreases as oc r~ 2 along the jet, i.e. the jet is 
isothermal. The emitting conical section of the jet may be thought of as the 
unresolved superposition of a number of shocks traveling down the jet, and 
shock acceleration keeps the nonthermal energy constant. Such a jet emits 
a flat radio spectrum up to the frequency where the spectrum steepens by 
one power due to the energy losses, typically in the submm-infrared regime. 
The flux at this break frequency is dominated by the jet cone near its apex. 
The threshold frequency z/ t h for the production of pions in head-on collisions 
between photons and protons with Lorentz factor 7 p is given by 



showing that the region of maximum surface brightness temperature is most 
important for the cooling of protons. Owing to electromagnetic cascading the 
electromagnetic power injected into the jet plasma by the cooling baryons is 
redistributed smoothly over the X-ray and gamma ray bands. Therefore, the 
emission component with the largest energy flux dominates the entire high- 
frequency spectrum and proton-initiated cascade spectra have been computed 
only for this zone. An integration along the jet axis (which is necessary to 
obtain the flat radio spectrum) would only lead to marginal corrections in the 




(6) 
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gamma ray regime. The geometry is depicted in Fig. 2 together with a sketch 
of the emission components from the various scales. Due to the small size of 
the region dominating the gamma ray energy flux compared with the total vol- 
ume of the conical jet, traveling shocks would make the gamma ray emission 
much more susceptible to flux variability than the radio- infrared emission in 
qualitative agreement with the observations. As long as the perturbation time 
scale is larger than the proton acceleration time scale, the emission spectra 
can still be computed as quasi-stationary spectra, and it is an open theoreti- 
cal challenge to solve for the spectra in the general time- dependent case. An 
interesting time-dependent solution exists for a simplified version of Compton- 
scattering dominated proton-initiated cascades (22). Using the Doppler factor 
5 to convert between comoving and observer's frame, the proton acceleration 
time scale assuming Bohm diffusion is given by 

in comfortable agreement with some of the observed TeV flux variation time 
scales. The comoving frame gyro-radius of the protons 

r ~ 10 16 ( JL) 1 ( -\) cm (8) 
g V30 GJ V10 11 / V ; 

satisfies the constraint r g < rj if the jet radius is the lightcylinder radius of an 
MHD jet from a black hole of mass M and Schwarzschild radius r s in which 
case 7j ~ 100r s — 3 x 1O 16 (M/1O 9 M ) cm. For time scales longer than the 
proton acceleration time scale, correlated flux variations from hard X-rays to 
TeV gamma rays are expected as a typical phenomenon (within the model 
assumptions). If the hard X-ray emission or the gamma ray emission in the 
EGRET band is dominated by emission from the accelerated electrons, the 
TeV variability can be different. This case would argue for a higher ratio of the 
photon to magnetic energy density u^/ub ~ 1 than in the generally assumed 
case u-y/uB < 1 in the proton blazar model. Note also that variability on 
shorter time scales must be expected under realistic conditions, e.g. due to the 
presence of inhomogeneities along the jet. These variations can occur on time 
scales down to the cooling time scale of the cascade electrons, and require the 
solution of the time-dependent cascade equations. The short-term variability 
behavior of the cascades is very complicated owing to the pair production 
threshold varying with gamma ray energy. 

The TeV detections of nearby blazars confirm the prediction of the proton 
blazar model that the relativistic jets should be optically thin to gamma ray 
emission below ~ 1 TeV. This means that the conversion of injected gamma 
rays into pairs and vice versa becomes unimportant below TeV (unsaturated 
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v b v c log v <\l h 

Fig. 2. Conical geometry for the relativistic jet with Lorentz factor 7j and opening 
angle <I> ~ V^Tj assumed in the proton blazar model. Note that the total pro- 
ton-initiated cascade emission from the cone is dominated by the emission from the 
region of highest surface brightness temperature in the submm-infrared regime at 
r = rb and is therefore more susceptible to variations due to traveling shocks (at 
position r s with velocity /3 S ). The sketch of the spectral energy distribution indicates 
the flux contributions from various scales r in the jet. An additional effect increas- 
ing the gamma ray flux from the region of highest surface brightness temperature is 
indicated by the parameter £ which is proportional to the proton maximum energy 
and which is expected to decrease with r due to the nonlinear development of the 
shock structure. 

cascades). Emission above ~ 1 TeV is expected to be optically thick, but 
nevertheless of approximate power law shape. Synchrotron cooling is assumed 
to be the dominant process replenishing the gamma rays from the pairs. For 
the entire parameter space, synchrotron emission remains non-relativistic so 
that the characteristic synchrotron photon energy is much less than that of the 
radiating electrons assuring rapid convergence of the cascade equations (which 
can be brought into the form of a Volterra integral equation of the second 
kind) by Banach's fixed point theorem. The shape of the multi-TeV spectrum 
follows from very simple considerations. It is assumed that the protons have 
a differential distribution dN/dE oc E~ Sp . The electron distribution has the 
same slope in the optically thin range, but in the energy range responsible 
for producing the target photons for the protons their spectrum is steeper by 
one power s e = s p + 1 owing to energy losses. The electrons thus produce a 
synchrotron flux density spectrum S u oc v~ a with spectral index a = s p /2. 
In the original papers (7) it was assumed that s p = 2 and correspondingly 
a = 1 which is the non-relativistic result from 1st order Fermi acceleration 
at strong shocks in the test-particle approximation. A flatter value s p = 3/2 
corresponding to a = 3/4 may be more appropriate for strong shocks in the 
general case (23). The cascades are initiated by gamma rays from the decay of 
the neutral pions at ultrahigh energies. The slope of the differential injection 
spectrum of the gamma rays is given by a, but steepens due to pair creation 
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on the synchrotron target characterized by the optical depth t(E) oc E a . 
The steepening can be described by the energy-dependent escape probability 
-Pesc = 1 — exp[— t(E)]/t(E) — > t(E)^ 1 oc E~ a for t>1, i.e. a steepening of 
the injection spectrum by a. Hence, the stationary injection spectrum has the 
slope s 7j i = 2a. The next step in the cascade development involves the creation 
of pairs which have a stationary distribution with the same slope 2a producing 
the second generation of gamma rays. Since generally these gamma rays still 
lie in the optically thick energy range, the stationary gamma ray distribution 
has the slope s 7)2 = 2a + 0.5. The same is true for the new generation of 
pairs produced by these gamma rays, and their synchrotron gamma rays are 
mostly emitted at optically thin energies where their spectrum has the slope 
s 7 ,3 = a + 0.75 (in the optically thick range the slope is 2a + 0.75). Thus 
the predicted multi-TeV slope is bracketed by 2a + 0.5 and 2a + 0.75. For 
a — 1, the corresponding range is 2.5 — 2.75 and for a = 0.75 it is 2.0 — 
2.25 which is in reasonable agreement with the observations of Mrk 421 and 
Mrk 501 (11) corrected for the expected intergalactic gamma ray attenuation 
(12). The model spectrum fitted to lower frequency data as published prior 
to the 1997 HEGRA multi-TeV observations agrees remarkably well with the 
measurements (10). The shape of the multi-TeV spectrum is not sensitive 
to changes in the maximum energy and can remain constant under large- 
amplitude changes of the flux associated with changes in the maximum energy. 



4 Neutrino and cosmic ray predictions 



The photo-production of charged pions leads to the emission of neutrons and 
neutrinos. Neutrons associated with the production of 7r + have no efficient cou- 
pling with the magnetized plasma in the jet and therefore escape ballistically. 
The neutrons decay to protons after a propagation length l n = (7 n /10 11 ) Mpc, 
and such extragalactic cosmic rays suffer energy losses traversing the mi- 
crowave background (25). At an observed energy of 10 19 eV, the energy- 
loss distance is A p ~ 1 Gpc owing to pair production. This distance cor- 
responds to a redshift z p determined by A p = (c/H ) Jq p dz/[{l + z)Ejp(z)] 

where E 3P (z) = + zf + Q R (1 + z) 2 + fi A ] 5 with Q + tt R + Q A = 1. 
Almost independent on cosmology, the resulting value for z p is given by 
z p — ^50/(6 — h 50 ) — 0.2/i 50 where h 50 = H o /50 km s -1 Mpc -1 . Therefore, 
when computing the contribution of extragalactic sources to the observed 
cosmic ray flux above 10 19 eV, only sources with z < z p must be consid- 
ered. Assuming further that extragalactic sources of cosmic rays and neutri- 
nos are homogeneously distributed with a monochromatic luminosity density 
^(z) oc (1 + z) 3+k where k ~ 3 for AGN (24), their contribution to the energy 
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density of a present-day diffuse isotropic background is given by 

,(0) = / + ,) '-rfz = ^/ L > (9) 

v ' y ' 



where z m = 2 denotes the redshift of maximum luminosity density. The factor 
(1 + z)~ A accounts for the expansion of space and the redshift of energy. For a 
simple analytical estimate of the effect of energy losses on the proton energy 
density at 10 19 eV, we collect only protons from sources out to the horizon 
redshift z p ~ 0.2 for 10 19 eV protons, whereas neutrinos are collected from 
sources out to the redshift of their maximum luminosity density z m . This yields 
the energy density ratio for neutrinos at an observed energy of ~ 5 10 17 eV 
and protons at 10 19 eV 

MO) _ Uo Zin (l + z) k - 2 /Ej P (z)dz 

«p(0) C(l + z)^/E JP (z)dz ~ 1 ) 



using £ ~ 0.3 from decay and interaction kinematics, and considering an open 
Universe with Ejp(z) = (I+2) and a closed one with Ejp(z) = (1+z) 3 / 2 . Fig. 3 
shows the exact propagated proton and neutrino spectra for Q = 1 from a full 
Monte-Carlo simulation employing the matrix doubling method of Protheroe 
& Johnson (26). The assumed neutron spectrum was dN n /dE n oc E~ x (cor- 
responding to a = 1 in the previous section) up to 10 18 eV and dN n /dE n oc 
E~ 2 exp(—E n /E cnt ) above. The steepening reflects the fact that the maximum 
energy may vary from source to source. The muon neutrino spectra follow the 
same shape, but they are shifted according to a simplified treatment of pion de- 
cay and production kinematics. A more accurate treatment yields small correc- 
tions (29). The neutron spectrum was normalized to the cosmic ray data yield- 
ing a neutrino spectrum consistent with model A from the original work (27). 
The associated gamma ray flux F 1 « 2F V ln[100] ~ 2 x 10~ 6 GeV cm -2 s _1 sr _1 
corresponds to a sizable fraction of the observed gamma ray background fluxQ- 
The neutrino flux is consistent with the bound given in ref. (28), although it is 
possible to have extragalactic neutrino sources of higher neutrino fluxes with- 
out violating the observed cosmic ray data as a bound (29). One could easily 
construct such models which produce the entire gamma ray background on 
the same rationale (29). 



2 A recent paper by Waxman and Bahcall (28) refers to the neutrino flux from 
model B in the original work which was given only to demonstrate that hadronic jets 
can not produce a diffuse gamma ray background with an MeV bump (as measured 
by Apollo and which is now known to be absent from a COMPTEL analysis) without 
over-producing cosmic rays at highest energies. 
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Fig. 3. Comparison of proton (solid line) and neutrino fluxes (dotted lines, from top 
to bottom iv, u^, and u e ) from the proton blazar model (Monte-Carlo computations 
and figure kindly provided by R.J. Protheroe). 



5 Neutrino oscillations and event rates 



The neutrino flux shown in Fig. 3 corresponds to a very low muon event rate 
even in a detector with an effective area of 1 km 2 (~ 1 event per year and 
per steradian above 100 TeV). This event rate is a very conservative estimate, 
since there must be additional neutrino production due to pp-interactions of 
escaping nucleons diffusing through the host galaxies and galaxy clusters. The 
neutrino flux could also be increased by increasing the number of extragalactic 
jet sources with proton maximum energies well below 10 19 eV (29). As a matter 
of fact, such a model ramification is required if one wants to explain the entire 
diffuse gamma ray background by hadronic photo-production sources. 

At this point the discovery of neutrino mass announced by the Super-Kamiokande 
collaboration (30) changes the situation in a major way. A deficit of atmo- 
spheric muon neutrinos was observed with Super-Kamiokande at large zenith 
angles with the most likely explanation being a full-amplitude oscillation of 
muon flavor eigenstates to tauon flavor eigenstates across the Earth at GeV 
energies. The transition probability -P(f M — > v T ) is a function of distance and 
energy L/E, i.e. 

P(v u -> u T ) « sin 2 20 sin 2 11 



where the mixing angle 9 specifies the mixing amplitude sin 2 (29) > 0.82. The 
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missing piece of information about the neutrino oscillations is the appearance 
of tau leptons to which Super-Kamiokande is not sensitive. A long-baseline ex- 
periment using muon neutrinos from a laboratory beam is extremely difficult if 
not impossible with existing laboratories, since one must establish a very large 
distance and a high energy to reach the tau mass shell at m T = 1.784 GeV. 
The beam luminosity decreases rapidly with energy and distance which poses 
an irreducible problem. However, an astrophysical beam of muon neutrinos, 
such as the one proposed in this paper and for which the multi-TeV obser- 
vations give us somewhat more confidence that they really exist, is ideally 
suited for this type of experiment, since both L and E obtain "astronomically 
large" values. Hence it follows that very likely astronomical high-energy neu- 
trino sources such as extragalactic radio jets will help to solve a major puzzle 
in elementary particle physics. 

There is another effect associated with neutrino oscillations which alleviates 
the problem with the low muon event rate of the predicted flux. To see this, 
one must realize that the solid angle for the detection of high-energy muon 
neutrinos becomes very narrow at energies in excess of ~ 100 TeV, since the 
Earth becomes optically thick to muon neutrinos above this energy (the weak 
interaction cross section depends on energy). This is called the Earth shad- 
owing effect (31). Below 100 TeV the atmospheric background of neutrinos is 
too strong for the discovery of the rare events due to extragalactic neutrinos, 
unless the angular resolution of the neutrino telescopes is very good. Actu- 
ally, the atmospheric background is comparable to the predicted astronomical 
background at around 100 TeV. Therefore one is confronted with the problem 
of a too small solid angle for events above 100 TeV and a too low rate below 
100 TeV. If the neutrino oscillation hypothesis is correct, the extragalactic 
muon neutrino beam is fully mixed with tau neutrinos at Earth. Tau leptons 
produced by charged-current interactions inside the Earth decay before further 
interacting. Since the Earth is also opaque to tau neutrinos above ~ 100 TeV, 
all the tau neutrinos entering on the one side with energies above 100 TeV 
emerge on the other side with energies of around 100 TeV (32) obliterating the 
Earth-shadowing effect. In the proton blazar model with a — 1, the number 
density of neutrinos above 100 TeV remains constant up to roughly 10 5 TeV 
implying an increase in the contained event rate by a factor 

1 °f dx 



fr = ^ J — = 111 L 10 J " 3 ' 5 

10 2 



where the factor \ is for the mixing between muon and tauon neutrinos (| 
would be appropriate for further mixing with electron neutrinos). The neu- 
trinos from extragalactic radio sources can therefore be expected to produce 
more events than the atmosphere at around 100 TeV. Some of these events are 
not muon events, but direct tauon events. Owing to tauon decay, the tauon 
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tracks are very short l T = 50(i? T /100 TeV) cm and the main signature is the 
electromagnetic cascade from tau decay. Horizontal events can be of higher en- 
ergy, producing the famous double-bang events with the first bang indicating 
the charged-current tauon production event and the second its decay. Thus, 
the glass is half-full. 



6 Discussion and summary 

The multi-TeV spectra from nearby blazars prediced on the basis of the proton 
blazar model are in accord with the observations if the effect of pair atten- 
uation due to the extragalactic infrared background is taken into account. 
Variability patterns are similar to those in synchrotron- self- Compton models, 
but more complex, since variations of the target photon flux fold into the 
cascade development in a non-linear manner. Short time scale variability is 
likely to reflect the passage of shocks through inhomogeneities and correspond 
to cooling time scale variations. If the cosmic ray flux emitted by hadronic 
accelerators is enough to explain the observed cosmic rays above 10 19 eV, 
the associated gamma ray power from these sources is enough to produce at 
least a sizable fraction of the observed extragalactic gamma ray background. 
The gamma ray power is larger than that in cosmic rays, since the cosmic 
rays lose energy traversing the low-energy background radiation fields and 
most sources have high redshifts. Strong evolution of their luminosity density 
would rule out GRBs as possible sources of the highest energy cosmic rays, 
since their cumulative gamma ray flux is far below the extragalactic gamma 
ray flux. Although the expected muon event rate in neutrino telescopes is low, 
neutrino oscillations lead to tau cascades canceling the Earth shadowing effect 
thereby increasing the detection probability. 



References 

[1] G.R. Burbidge, Astrophys. J. 124 (1956) 416. 
[2] G.C. Perola, Astron. & Astrophys. 3 (1969) 481. 

[3] R.D. Blandford, R.L. Znajek, Mon. Not. Roy. Astr. Soc. 179 (1977) 
433. 

[4] J.F.C. Wardle, et al., Nature 395 (1998) 457. 

[5] P.L. Biermann, P.A. Strittmatter, Astrophys. J. 322 (1987) 643. 

[6] M. Sikora, et al, Astrophys. J. 320 (1987) L81. 

[7] K. Mannheim, W. Kriills, P.L. Biermann, Astron. & Astrophys. 

222(1991) 222, K. Mannheim, P.L. Biermann, Astro. & Astrophys. 

333(1992) L21, K. Mannheim, Astron. & Astrophys. 269 (1993) 67. 
[8] P. Meszaros, M.J. Rees, Astrophys. J. 405 (1993) 278. 



13 



[9] R. McCray, Astrophys. J. 156 (1969) 329. 
[10] K. Mannheim, Science 279 (1998) 684. 
[11] A. Konolpenko, et al. (HEGRA Collaboration), this volume 
[12] O.C. de Jager, E. Dwek, Astropart. Phys. (1999), submitted; F.W. 

Stecker, O.C. de Jager, Astron. & Astrophys. 334 (1998) L85. 
[13] M.J. Rees, J. Silk, Astron. & Astrophys. (1998) 331, LI. 
[14] D.B. Sanders, et al., Astrophys. J. 347 (1989) 29. 
[15] D.E. Gruber, The X-ray Background, eds. X. Barcons and A.C. Fabian, 

Cambridge U.R (1992) p. 44. 
[16] S. Rawlings, S. Saunders, Nature 349 (1991) 138. 

[17] P. Sreekumar, F.W. Stecker, S.C. Kappadath, in: Proceedings of the 
Fourth Compton Symposium, eds. Charles D. Dermer, Mark S. Strick- 
man, and James D. Kurfess, Williamsburg, VA, April 1997: AIP Con- 
ference Proceedings 410 (1998) p. 344. 

[18] P. Madau, E.S. Phinney, Astrophys. J. 456 (1996) 124. 

[19] M.H. Salamon, F.W. Stecker, Astrophys. J. 493 (1998) 547. 

[20] H. Bohringer, W. Voges, A.C. Fabian, A.C. Edge, D.M. Neumann, 
Month. Not. Roy. Astr. Soc. 264 (1993) L25. 

[21] R.D. Blandford, A. Konigl, Astrophys. J. 232 (1979) 34. 

[22] A. Mastichiadis, J. Kirk, Astron. Astrophys. 295 (1995) 613. 

[23] M.A. Malkov, Astrophys. J. 511 (1999) L53. 

[24] B. J. Boyle and R. J. Terlevich, Mon. Not. Roy. Astro. Soc. 293 (1998) 
L49. 

[25] J. P. Rachen and P. L. Biermann, Astron. & Astrophys. 272 (1993) 161. 

[26] R. J. Protheroe and P. A. Johnson, Astropart. Phys. 4 (1996) 253. 

[27] K. Mannheim, Astropart. Phys. 3 (1995) 295. 

[28] E. Waxman, J. Bahcall, Phys. Rev. D (1999) accepted. 

[29] K. Mannheim, R.J. Protheroe, J.P. Rachen, Phys. Rev. D. (1999) sub- 
mitted; F.W. Stecker, M.H. Salamon; Space Set. Rev. 75 (1996) 341. 

[30] Y. Fukuda, et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 
81 (1998) 1562. 

[31] R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, Astropart. Phys. 5 
(1996) 81. 

[32] F. Halzen, D. Saltzberg, Phys. Rev. D. 81 (1998) 4305. 



14 



